A&A 546, A1 (2012) 

DOI: 10.1051/0004-6361/201219718 
© ESO 2012 


Astronomy 

Astrophysics 


Chandra grating spectroscopy of three hot white dwarfs 

J. Adamczak 1 ’ 2 , K. Werner 1 , T. Rauch 1 , S. Schuh 3 , J. J. Drake 4 , and J. W. Kruk 5 


1 Institute for Astronomy and Astrophysics, Kepler Center for Astro and Particle Physics, Eberhard Karls University, Sand 1, 
72076 Tubingen, Gennany 

e-mail: adamczakOastro .as.utexas.edu 

2 McDonald Observatory, The University of Texas, Austin, Texas, 78712, USA 

3 Institute for Astrophysics, Georg August University, Friedrich-Hund-Platz 1, 37077 Gottingen, Germany 

4 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 

5 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 

Received 30 May 2012 / Accepted 8 August 2012 


ABSTRACT 

Context. High-resolution soft X-ray spectroscopic observations of single hot white dwarfs are scarce. With the Chandra Low-Energy 
Transmission Grating, we have observed two white dwarfs, one is of spectral type DA (LB 1919) and the other is a non-DA of spectral 
type PG 1159 (PG 1520+525). The spectra of both stars are analyzed, together with an archival Chandra spectrum of another DA 
white dwarf (GD 246). 

Aims. The soft X-ray spectra of the two DA white dwarfs are investigated in order to study the effect of gravitational settling and 
radiative levitation of metals in their photospheres. LB 1919 is of interest because it has a significantly lower metallicity than DAs 
with otherwise similar atmospheric parameters. GD 246 is the only white dwarf known that shows identifiable individual iron lines in 
the soft X-ray range. For the PG 1159 star, a precise effective temperature determination is performed in order to confine the position 
of the blue edge of the GW Vir instability region in the HRD. 

Methods. The Chandra spectra are analyzed with chemically homogeneous as well as stratified NLTE model atmospheres that assume 
equilibrium between gravitational settling and radiative acceleration of chemical elements. Archival EUV and UV spectra obtained 
with EUVE, FUSE, and HST are utilized to support the analysis. 

Results. No metals could be identified in LB 1919. All observations are compatible with a pure hydrogen atmosphere. This is in 
stark contrast to the vast majority of hot DA white dwarfs that exhibit light and heavy metals and to the stratified models that 
predict significant metal abundances in the atmosphere. For GD 246 we find that neither stratified nor homogeneous models can fit 
the Chandra spectrum. The Chandra spectrum of PG 1520+525 constrains the effective temperature to T e g = 150000 + 10 000 K. 
Therefore, this nonpulsating star together with the pulsating prototype of the GW Vir class (PG 1 159-035) defines the location of the 
blue edge of the GW Vir instability region. The result is in accordance with predictions from nonadiabatic stellar pulsation models. 
Such models are therefore reliable tools to investigate the interior structure of GW Vir variables. 

Conclusions. Our soft X-ray study reveals that the understanding of metal abundances in hot DA white dwarf atmospheres is still 
incomplete. On the other hand, model atmospheres of hydrogen-deficient PG 1159-type stars are reliable and reproduce well the 
observed spectra from soft X-ray to optical wavelengths. 

Key words, white dwarfs - stars: abundances - stars: atmospheres - stars: individual: LB 1919 - stars: individual: GD246 - 
stars: individual: PG 1520+525 


1. Introduction 

An unexpectedly small number of white dwarfs (WDs) were 
detected in the ROSAT PSPC X-ray all sky survey (Fleming 
et al. 1996). Most of them (161) are of spectral type DA (pure- 
hydrogen optical spectra), with an additional three DOs (helium- 
dominated), three DAOs (mixed H/He optical spectra), and eight 
PG1159 stars (He-C-0 dominated). It was realized that the 
atmospheric opacity of radiatively levitated metals effectively 
blocks the outward leakage of X-rays in the vast majority of DAs 
with r e ff > 40 000 K (Barstow et al. 1997). Furthermore, the in- 
terstellar medium proved denser than expected in many lines of 
sight (see, e.g., Welsh et al. 1999), thus affecting cooler WDs. 

High-resolution soft X-ray spectroscopy of hot WDs enables 
the identification of chemical species in their photospheres that 
cannot be detected in other wavelength ranges. However, such 
observations are rather scarce and only became feasible with 
the advent of the Chandra observatory with its Low-Energy 
Transmission Grating (LETG). 


Until recently, only five WDs were observed with Chandra 
LETG. H 1504+65 is an extremely hot ( T e g = 200000 K), pecu- 
liar PG 1159 star, obviously a naked C-0 or O-Ne WD, whose 
soft X-ray spectrum is characterized by highly ionized O, Ne, 
and Mg lines (Werner et al. 2004a). On the other hand, GD 246 
is a DA, and its Chandra spectrum allowed the first unambigu- 
ous identification of lines from highly ionized iron (Vennes & 
Dupuis 2002). The Chandra spectra of the two DAs HZ 43 and 
Sirius B are perfectly matched by pure-hydrogen atmospheres 
and are used as soft-X-ray calibration targets (Pease et al. 2000; 
Vennes & Dupuis 2002; Beuermann et al. 2006). A Chandra ob- 
servation of the hot DO KPD 0005+5 106 (T e f f = 200 000 K) was 
used to prove that the soft X-ray emission is of photospheric and 
not coronal origin; however, the signal-to-noise ratio (S/N) is in- 
sufficient to identify individual spectral lines (Drake & Werner 
2005). 

Few other WDs are bright enough in soft X-rays to 
obtain useful Chandra spectra. We have observed two of 
them (LB 1919 and PG 1520+525) and including archival data 
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of GD246 present a spectral analysis in this paper. The 
DA LB 1919 is of interest because previous Extreme Ultraviolet 
Explorer (EUVE) spectroscopy showed an unexpectedly high 
soft X-ray flux that reveals an unexplained low metal abundance 
(Wolff et al. 1998), in contrast to expectations from radiative- 
levitation theory. We performed X-ray spectroscopy in order to 
identify individual elements that could hint at the origin of the 
metal deficiency. 

We observed PG 1520+525 in order to constrain its effec- 
tive temperature. The motivation lies in the fact that this non- 
pulsating PG1159 star, together with the pulsating prototype 
PG 1159-035, defines the blue edge of the GWVir instability 
region in the Hertzsprung-Russell diagram (HRD), see e.g. Jahn 
et al. (2007). The homogeneous atmosphere of PG 1520+525 
provides a convenient comparison for the investigated DA white 
dwarfs because it is not subject to the uncertainties in the treat- 
ment of radiative levitation and stratification of atmospheric 
composition. Thus, PG 1520+525 provides a test case for the 
model atmospheres employed in this study. A failure to achieve 
a consistent fit to the UV/optical and X-ray spectra of the ho- 
mogeneous PG1159 photosphere would indicate serious model 
deficiencies in addition to problems that would be encountered 
with radiative levitation physics in the case of hot DA WDs. 

This paper is organized as follows. We first describe in more 
detail the motivation of our soft X-ray analyses (Sect. 2). Then 
we specify our model atmosphere calculations and the atomic 
data used (Sect. 3). In Sect. 4, we report on the observations uti- 
lized in our analysis. In Sects. 5 to 7 we delineate the analysis 
procedure of our three program stars one at a time, and in Sect. 8 
we conclude. 


2. Soft X-ray emission from white dwarfs 

A concise presentation of past soft X-ray observations of WDs 
and their interpretation can be found, e.g., in Werner (2008). 

2. 1. DA white dwarfs 

DA WDs possess almost pure hydrogen atmospheres. The high 
gravity in these objects results in a chemical stratification with 
the lightest element, hydrogen, floating atop. In hot DAs {T t ff > 
20 000 K), hydrogen is almost completely ionized and the opac- 
ity in the atmosphere is strongly reduced. Soft X-ray radia- 
tion can emerge from deep, hot, photospheric layers. However, 
this radiation may be blocked by the opacity of heavy ele- 
ments that can be kept in the atmospheres by radiative levita- 
tion (Chayer et al. 1995a, b). The increasing efficiency of this 
mechanism with T e ff results in only a few X-ray detected DAs 
with J eff > 60 000 K. 

Since the spectral resolution of EUVE is too low to identify 
lines of individual species, the metal abundances of the hitherto 
investigated DAs were determined relative to the well-studied 
G 191-B2B (Wolff et al. 1998) using chemically homogeneous 
model atmospheres. The scaling of the metallicity in these mod- 
els relative to that of G 191-B2B (the so-called metallicity in- 
dex) results in satisfying spectral energy distribution (SED) fits 
for many, but not all, DAs. 

To make further progress, non-LTE model atmospheres were 
developed that calculate the abundances of the elements at each 
atmospheric depth in a self-consistent way, assuming equilib- 
rium between gravitational downward pull and radiative upward 
acceleration (Dreizler 1999). The EUVE SEDs of most DAs can 
be reproduced well by these models (Schuh et al. 2002). Some 


stars show a higher metallicity than predicted. This can be ex- 
plained, for instance, by ongoing accretion of circumstellar or 
interstellar matter. Other stars have a lower metallicity than pre- 
dicted by these models, and the reason for this is not known. A 
drastic example is HZ43A ( T e fs = 51 000 K, Beuermann et al. 
2006) whose SED is perfectly matched by a pure H model atmo- 
sphere. Another low-metallicity DA is the even hotter LB 1919 
(T e s = 56 000 K, this paper), which is subjected to analysis here. 

The exact knowledge of element abundances in the atmo- 
spheres of the metal-poor DAs might give hints as to an expla- 
nation for their metal deficiency. Many spectral lines of the high 
ionization stages of heavy elements like iron and nickel lie in 
the soft X-ray wavelength range. Spectroscopy with Chandra 
enables us to investigate these objects more closely. Our work 
represents the first detailed analysis of high-resolution soft X-ray 
spectra of DAs. Besides LB 1919, we study the more metal-rich 
DA GD 246 as a reference object. 

2.2. PG 1 159 stars 

As mentioned in the introduction, only very few non-DAs (DO 
and PG1159 stars) were detected in the soft X-ray band be- 
cause of the additional opacity of helium and enriched met- 
als. PG 1159 stars are hot, H-deficient (pre-) WDs with atmo- 
spheres mainly composed of He, C, and O. It is thought that 
their H-envelope was consumed during a late He-shell flash 
(Werner & Herwig 2006, and references therein). Among the 
PG 1159 stars, PG 1520+525 (r eff = 150000 K) is one of the 
brightest soft X-ray sources and hence an interesting target to 
study the spectral characteristics of this class in that wavelength 
region. The surface abundances of these hot stars are not affected 
by gravitational and radiative acceleration because of a weak 
radiation-driven wind (Unglaub & Bues 1997). Consequently, 
their atmospheres and envelopes can be assumed to be chemi- 
cally homogeneous. 

Our primary aim for the Chandra observation of the nonvari- 
able PG 1520+525 is to constrain its effective temperature. Then 
we can compare its position in the log T e ff - log g diagram with 
that of the prototype of the GW Vir pulsators, PG 1159-035. 
The blue edge of the instability region is confined by these stars 
(Werner et al. 1996). Based on nonadiabatic computations, the 
most advanced pulsation models for GW Vir stars, are those pre- 
sented by Quirion et al. (2007). Accordingly, the exact location 
of the blue edge depends on the envelope composition, primar- 
ily on the C and O abundance and, to a smaller extent, on the 
metallicity. Therefore, the concept of a blue edge is necessarily 
“fuzzy”. It was shown by spectroscopic analyses that both stars 
have (within error limits) the same abundances of C, O, and Fe 
(see below) and thus define the blue edge for that particular com- 
position. We can therefore prove or disprove the prediction of 
corresponding pulsation models as to the location of the edge. 
This represents a strong test for these models and their potential 
to derive the interior structure of PG 1159 stars by asteroseis- 
rnology methods. 

3. Model atmospheres and atomic data 

For the spectral analysis, chemically homogeneous NLTE 
model atmospheres were computed with the Tubingen Model 
Atmosphere Package (TMAP 1 , Werner et al. 2003), and chem- 
ically stratified models with a special variant, the diffu- 
sion/levitation NGRT code (Dreizler 1999). In comparison to 

1 http : //astro . uni - tuebingen . de/~TMAP 
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earlier work with NGRT, refinements to the code were applied 
for more realistic physics. Moreover, new model atoms for Fe, 
Ni, and Ge were constructed. These improvements are described 
in the following. 

3.1. Improvements to the diffusion/levitation code NGRT 

In the NGRT code, the element abundance pattern at any depth 
point in the atmosphere is calculated self-consistently by assum- 
ing equilibrium between gravitational and radiative forces. As a 
result, the atmosphere is no longer chemically homogeneous but 
vertically stratified. 

3.1.1. Bound-free transitions 

The primary assumption for the calculation of the radiative force 
of a bound-free (bf) transition in NGRT is that the photon mo- 
mentum is transferred completely to the ion, keeping the elec- 
tron out of consideration. As a consequence, the calculated ra- 
diative acceleration of the remaining ion might be overestimated. 
To check this assumption quantitatively, a correction according 
to Gonzalez et al. (1995) was implemented and tested. 

The radiative force on a particle of element A, with ionic 
charge i for a bound-bound (bb) or bf transition j, caused by 
photons in the frequency range (v, v + dv) is 


x 

X 

05 

o 


X 

x 

05 

o 

< 





Fig. 1 . Top panels: element abundances over depth following from 
two different treatments of radiative acceleration by photoionization 
(r e ff = 54 000 K, logg = 8.2). Left : the whole photon momentum is 
transferred to the ion. Right: only a fraction is transferred, according to 
Eq. (3). Bottom panel: difference. The peaks correspond to a shift of the 
steep regions of the abundance pattern by one depth point in the model 
atmosphere. 


3.1.2. Bound-bound transitions 


F'J dv = -j^-dv = <Tij(v)— dv, 


dp 
d t 


( 1 ) 


with the net momentum d p v , transported by the radiation flux T y . 
As mentioned above, in the case of a bf transition, the momen- 
tum of the photon causing an ionization of ion A' is not com- 
pletely transferred to ion A' +1 . A part of it is taken away by 
the ejected electron. Thus, a correction factor f oa has to be in- 
troduced, describing the remaining fraction of the momentum 
transferred to the ion A' + 1 . 

The correction factor can be expressed in terms of the fre- 
quency v of the photon, the threshold frequency vo necessary to 
eject the electron and a factor a\. 


4 (v - v 0 ) 

/ion — 1 0 a \‘ 

3 v 


(2) 


For testing purposes, the quantum calculation of Sommerfeld 
(1939) was taken and applied to the code. With his value of 
fli = 6/5, /ion becomes 


8 hv - x 
5 hv 


( 3 ) 


with the ionization energy threshold^. An interesting result of 
Eq. (3) was already stated by Michaud (1970). For certain fre- 
quencies the electron can be ejected with more momentum than 
brought in by the photon (/ e i ec = 1 -/i on > 1). As a consequence, 
the ion will be pushed back in the atmosphere by the photoion- 
ization. Thus, including a detailed bf absorption calculation can 
in principle lead to either higher or lower outward directed ra- 
diative forces on an ion. 

The comparison between abundances at each depth of the 
atmospheres calculated with old and new treatments of bf tran- 
sitions reveals, however, only minor and practically unimportant 
deviations (Fig. 1); there is no effect on the SED, in particular on 
the spectral lines, of the model. 


In the original NGRT code as used by Schuh et al. (2002), the 
redistribution of transferred photon momentum over the ion- 
ization stages A' and A i+1 by a bb transition in A' was treated 
such that the momentum was completely transferred to the next 
higher ion A ,+1 . This treatment is based on the assumption that 
an ionization following the radiative excitation of the ion takes 
place before the ion is deexcited by an inelastic electron colli- 
sion. Especially in very hot atmospheres, this need not be true. 
For a realistic treatment, the probabilities for ionization and col- 
lisional deexcitation would have to be evaluated. In a statisti- 
cal sense, fractions of the momentum are transferred to the next 
higher ion A 1+1 and to the originally excited ion A'. 

As a rough approach to test the maximum possible effect, 
the calculation of the bb transitions was modified such that the 
entire photon momentum is transferred to the lower ion A'. For 
example, we show the effect of the modified bb (and bf) pro- 
cedures on the population numbers of iron superlevels (intro- 
duced below) in Fig. 2. While the modified bf transitions leave 
the occupation numbers essentially unaltered, the modification 
of the bb transitions results in an obvious change, although only 
in the outer photospheric layers. The effect on the emergent flux 
is below 10% and thus considered insignificant, even under the 
extreme assumption of the fully inverted momentum distribu- 
tion. Therefore, a more complex treatment is unnecessary for 
our purposes. 

3.2. Model atoms 

For our model atmosphere calculations, we used two different 
sets of model atoms for the analysis of the PG 1 159 star on the 
one hand and for the two DAs on the other hand. The first set is 
identical to that described in detail by Werner et al. (2004a) in 
their Chandra analysis of H 1504+65. It comprises He, C, O, Ne, 
and Mg. The second set incorporates H, He, C, O, Si, P, S, and 
Ge, being represented by detailed model atoms taken from the 
Tubingen Model Atom Database 2 TMAD. It is augmented by Fe 

2 http : //astro . uni - tuebingen . de /~TMAD 
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Fig. 2. Effect of modified treatments of acceleration by line and photoionization transitions on occupation numbers of Fe v and Fe vi superlevels 
(r e ff = 56 000 K, logg = 7.9). The six panels show: 1. Original treatment of bb and bf transitions. 2. New treatment of bf transitions. 3. Difference 
between the first two. 4. New procedure for bb transitions. 5. bf and bb refinements are implemented. 6. Difference between the last two. The 
legend in panel 1 refers to the superlevel designations. 


and Ni which are treated in a statistical way using superlevels 
and superlines (Rauch & Deetjen 2003). Table 1 summarizes the 
set of data used for the DA analyses. 

The Ge model atom is new and presented in more detail 
below. We also use new Fe and Ni model atoms that comprise 
vastly extended line lists (Kurucz 2009) compared to all previ- 
ous analyses. Their characteristics are also described in detail 
below. 

3.2.1 . Germanium 

Germanium (Ge iv) was first identified in three hot DAs, among 
them GD 246, by Vennes et al. (2005). For GD 246, we perform 
the first NLTE analysis of Ge in a stellar atmosphere and in- 
vestigate its radiative levitation properties. The model atom was 
constructed using level energies from the NIST database. Due 
to the lack of data, all oscillator strengths fij are approximated 
by adopting values from the isoelectronic Civ, with the excep- 
tion of the two observed Geiv AA 1189.07, 1229. 84A lines, for 
which we use the fj from Chayer et al. (2005). Very recently, 
new fij for Ge iv lines were published (Nath Dutta & Majumder 
201 1). A single test calculation was performed in which we use 
the new data for all line transitions in the model atom. The effect 
on the observed two lines is insignificant. Photoionization rates 
were computed with hydrogen-like cross sections. Electron col- 
lisional excitation and ionization rates were evaluated with usual 
approximation formulae following van Regemorter (1962) and 
Seaton (1962), respectively. 

In order to check for the importance of NLTE effects, we per- 
formed another test calculation in which we enforced LTE popu- 
lations for Ge levels by artificially increasing electron collisional 
rates. The NLTE line profiles are significantly shallower than 
the LTE profiles (Fig. 4). As a consequence, the derived NLTE 
abundances will be roughly a factor of two higher compared to 
LTE analyses. 


3 http://physics.nist.gov/ 
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Fig. 3. Grotrian diagram of Ge iv. The resonance doublet is labeled. 



Fig. 4. Ge IV line profiles from LTE (dashed) and NLTE (solid) H+Ge 
model atmospheres (r eff = 57 000 K, log g = 7.6, Ge/H = 2.5 x 10 -9 ). 

3.2.2. Iron and nickel 

The Iron Opacity Interface (IrOnIc) constructs model atoms 
for iron-group elements as input for the stellar-atmosphere code 
(Rauch & Deetjen 2003). During the course of this work, new Fe 
and Ni line data were provided by Kurucz (2009) 4 . They include 

4 http : //kurucz . harvard . edu/atoms . html 
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Table 1. Statistics of model atoms used in the calculations for LB 1919 
and GD 246. 


Elem. 

Ion 

NLTE levels 


Lines 


H 

1 

10 


45 



II 

I 


- 


C 

II 

16 


37 



III 

58 


329 



IV 

54 


295 



V 

I 


- 


O 

II 

16 


26 



III 

21 


38 



IV 

18 


39 



V 

30 


95 



VI 

54 


291 



VII 

1 


- 


Si 

IP 

20 


48 



III 

17 


28 



IV 

16 


44 



V 

15 


20 



VI 

1 


- 


P 

III 

3 


0 



IV 

15 


9 



V 

18 


12 



VI 

1 


- 


S 

IIP 

3 


0 



IV 

17 


32 



V 

19 


32 



VI 

18 


48 



VII 

1 


- 


Ge b 

III 

14 


0 



IV 

8 


8 



V 

9 


0 



VI 

1 


- 


Fe 

III 

7 

(2467) 

27 

(537689) 


IV 

7 

(6389) 

27 

(3102371) 


V 

7 

(6728) 

25 

(3266247) 


VI 

8 

(5464) 

33 

(991935) 


VII 

9 

(2690) 

41 

(200455) 


VIII 

8 

(694) 

33 

(19587) 


IX 

1 


- 


Ni 

III 

5 

(3739) 

14 

(1033920) 


IV 

7 

(5939) 

27 

(2512561) 


V 

7 

(6165) 

25 

(2766664) 


VI 

7 

(10985) 

27 

(7408657) 


VII 

8 

(11866) 

33 

(4195381) 


VIII 

7 

(6469) 

27 

(1473122) 


IX 

1 


- 



Notes. Numbers in brackets denote levels and lines combined into su- 
perlevels and superlines. (o) LB 1919 models only. <4) GD246 models 
only. 

a much larger number of lines (about 20 times more), especially 
in the X-ray range. Figures 5 and 6 demonstrate the effect of 
these additional lines used in the atmosphere calculations. The 
new line lists were not incorporated in the bulk of the model cal- 
culations, however, the best-fitting models presented here were 
recalculated with the newest Kurucz atomic data sets available. 

4. Observations, reddening, interstellar H 
and He column densities 

We performed Chandra spectroscopy of LB 1919 and 
PG 1520+525 with the LETG/HRC-S. For GD246, we 
used the Chandra observation investigated by Vennes & Dupuis 
(2002). Other observations used in our work were retrieved 
from the MAST^ archive. All observations are listed in Table 2. 

5 http://archive.stsci.edu/ 



MA] 

Fig. 5. Comparison of H+Fe model spectra calculated with old and 
new Kurucz line lists (top left and right panels, respectively; T c g = 
57 000 K, log g = 7.90, Fe/H = 1.25 x 10 -6 ). Bottom panels', same 
models after convolution with the Chandra LETG instrument response 
(bin size 0.1 A). 



Fig. 6. Like Fig. 5, but here for nickel (Ni/H = 7.9 x 10 8 ). 


Table 2. Observation log of spectra analyzed in our study. 


Star 

Instrument (Id.) 

Observation 
start date (UT) 
mm-dd-yyyy 

^exp 

(ks) 

LB 1919 

EUVE 

04-16-1994 

166 


FUSE LWRS 

05-03-2001 

33 


IUE SWP 52808 

11-15-1994 

4.5 


Chandra LETG/HRC-S 

02-02-2006 

111 

GD 246 

EUVE 

08-08-1994 

14 


FUSE LWRS 

11-12-2000 

1.5 


FUSE MDRS 

07-14-2001 

24 


HST STIS/E140H 

11-20-1998 

2.4 


IUE SWP 40467 

12-27-1990 

0.33 


IUE LWP 19485 

12-27-1990 

0.84 


Chandra LETG/HRC-S 

01-14-2000 

40 

PG 1520+525 

Chandra LETG/HRC-S 

04-04-2006 

142 


In order to determine interstellar reddening, FUSE, IUE, and 
HST-STIS spectra (the latter marked with “S” in Fig. 7) as well 
as optical and IR magnitudes from the SIMBAD and NOMAD 
databases were investigated. Fluxes of our final models were 
scaled to the reddest available photometric magnitude. The red- 
dening law of Fitzpatrick (1999) was used with R v = 3.1. We 
determined E B ~v = 0.030 ± 0.005 and 0.0011 ± 0.0001 for 
LB 1919 and GD 246, respectively. These quantities will be used 


Al, page 5 of 13 


A&A 546, A1 (2012) 



Fig. 7. Determination of interstellar reddening of LB 1919 and GD 246 
from the observed flux distribution compared to models (thick red 
lines). For clarity, all spectra are convolved with 2 A (FWHM) 
Gaussians. 

below for distance determination. Accounting for absolute-flux 
uncertainties in the UV, more realistic errors are of the order 
±0.05. However, the error in the distance determination is domi- 
nated by the uncertainties in T e g and log g, and we neglect the er- 
rors for reddening. For PG 1520+525, Dreizler & Heber (1998) 
derived E B ~v = 0.0. 

Column densities for interstellar neutral hydrogen as well 
as neutral and ionized helium affect the analysis of EIJVE and 
Chandra spectra. 

LB 1919: Based on the values found in the EUVE analysis 
by Landenberger-Schuh (2005) (V(Hi) = 1.6 x 10 19 cm -2 , 
(V(Hei) = 0.04 V(Hi), V(Heii) = 0.04 JV(Hi)), the interstel- 
lar H and He column densities were varied in our fitting pro- 
cedure. Our values for two of the best-fitting models are de- 
picted in Fig. 12. For our present fit to the Chandra spectrum, 
we used V(Hi) = 1.6 x 10 19 cm” 2 , V(Hei) = 0.02 • V(Hi), 
V(Hen) = 0.04 ■ V(Hi). 

GD246: The Hi column density of 1.288 x 10 19 cm -2 was 
taken from Oliveira et al. (2003) and kept fixed. Our adopted 
values for the column densities of Hei (0.12 V(Hi)) and Hen 
(0.02 V(H i)) are close to the ones determined by Vennes et al. 
(1993) (V(Hei) = 1.05-1.25 x 10 18 cm” 2 , V(Hen) = 3.4-4.0X 
10 17 cm -2 ). 

Since these numbers depend on their assumed photospheric 
model, we adjusted the values slightly to better fit our mod- 
els. For our Chandra analysis presented here we chose 0.05 and 
0.026 MHi), respectively, which is close to the values used by 
Landenberger-Schuh (2005). 

PG 1520+525: Werner et al. (1996) used N( Hi) = 1.5 x 
10 2 °cm -2 for their fit of a T e g = 150 000 K model to the 
SED observed by EUVE. Dreizler & Heber (1998) determined 
(V(Hi) = 2.5 x 10 20 cm -2 from the Ly a line profile. For the fit 
to the Chandra spectrum, the Hi column density was treated 
as a free parameter. For our final model fits in this paper, we 
found 1.0 x 10 2 °cm -2 . The variation of the He column den- 
sities had no effect on the fit quality, and were kept fixed at 
V(Hei) = 0.05 V(Hi) and N( Hen) = 0.026- V(Hi). 

5. LB 1919 

5.1. Previous investigations 

LB 1919 (WD 1056+516) was analyzed by Vennes et al. (1997) 
using Balmer-line spectroscopy and pure-hydrogen model 



Fig. 8. Fit to the Lyman lines of LB 1919 (thin line). Overplotted are 
pure H models with T e s = 56 000 K, log g = 7.9 (thick line), and with 
= 69 000 K, log g =7.9 (dashed line). 


atmospheres. They obtained Teff = 68 640 K and log g = 
8.08 (cgs units). Finley et al. (1997) excluded the star from 
their analysis. Their spectrum exhibited flat-bottomed Balmer- 
line profiles, and they suggested orbital or rotational veloci- 
ties of ~ 1000 km s -1 as a possible interpretation, although they 
may be artifacts (Koester, priv. comm.). More recently, Lajoie & 
Bergeron (2007) determined T e g = 67 022 K and log g = 7.99, 
and Gianninas et al. (20 1 1 ) found T e g = 68 5 1 0 K and log g = 
7.94 by pure-hydrogen model fits to the Balmer lines. 

Wolff et al. (1998) investigated a sample of 20 DAs from 
the EUVE archive. Introducing the above-mentioned metallicity 
index (Sect. 2.1), they were able to reproduce the spectra of most 
of them. For LB 1919, they determined a metallicity at least ten 
times smaller compared to that of G 191-B2B, adopting T e g and 
log g values from Vennes et al. (1997). 

Based on these results, Landenberger-Schuh (2005) per- 
formed an analysis of EUVE spectra of DAs using new chemi- 
cally stratified NLTE model atmospheres with equilibrium abun- 
dances. In the case of LB 1919, the use of these new models 
resulted, however, in a worse fit compared to the analysis of 
Wolff et al. (1998) with homogeneous models because of too 
much levitation of metals in the diffusion models. The reason 
for this overprediction remained unexplained. 

5.2. FUSE 

A grid of pure-hydrogen model atmospheres was used to fit 
the Lyman lines in the FUSE spectral range. We found T ef j = 

56 000 ± 2000 K and log g = 7.9 ± 0.3. Somewhat hotter mod- 
els with higher gravity (e.g., T e g = 62 000 K, log g = 8.2) 
fit almost equally well, but are excluded by the Chandra spec- 
trum (Sect. 5.4). The temperature is significantly lower than 
the Balmer-line results reported above. It is remarkable that 
Lajoie & Bergeron (2007) arrived at a similarly low value of 

57 701 ± 1 8 400 K from the overall UV (IUE)/optical flux dis- 
tribution, and at an even lower value of 39 045 ± 18000K by 
fitting the Ly a line in the IUE spectrum. Figure 8 shows the 
comparison between a hotter model and a better-fitting, cooler 
model together with the FUSE Lyman lines. Adding metals to 
the 56 000 K model with abundances at the upper limit derived 
below has no significant influence on the Lyman-line profiles. 

Upper limits to metal abundances were determined with ho- 
mogeneous model atmospheres. Their accuracy is estimated to 
be about 0.3 dex following from the error ranges in T e g and log g. 
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Fig. 9. Homogeneous (thick line) and stratified (dashed line) model 
spectra near the Chi multiplet at A 1175 A in LB 1919. The homoge- 
neous model has C/H = 1.0 x 10 -7 . 
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Fig. 12. Best-fit homogeneous (solid) and stratified (dashed) models for 
the EUVE spectrum of LB 1919 (thin line). The homogeneous model 
has metal abundances set to the upper limits given in Table 4; neither 
model has oxygen. 


At T e j = 56 000 K, the stratified models predict no detectable 
lines only if the surface gravity is log g = 8.2 or higher. 


Fig. 10. Homogeneous (thick line) and stratified (dashed line) model 
profiles of P v and Siiv compared to LB 1919. The homogeneous model 
has Si/H = 5.0 x 10" 9 and P/H = 1.0 x KT 10 . 
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Fig. 11. Homogeneous (thick line) and stratified (dashed line) model 
profiles of S iv and Si iv compared to LB 1919. The homogeneous model 
has S/H = 1.75 x KT 8 and Si/H = 5.0 x 10“ 9 . 
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5.3. EUVE 

The stratified models calculated by Landenberger-Schuh (2005) 
had a flux too low to match the EUVE observation. The failure 
to derive a good fit resulted from the excessive opacity caused 
by the large number of elements incorporated in the model at- 
mospheres. We tried to reproduce the EUVE observation with 
stratified models that include only those metals listed in Table 4, 
i.e., metals that have been detected in the FUSE spectra of com- 
parable WDs, like GD246. Our results show that the strongest 
influence on the EUV SED comes from the opacity of oxygen 
and that a satisfying fit can only be derived for a model atmo- 
sphere without this species and at high gravity (logg = 8.5). We 
also computed homogeneous models with the same elements and 
abundances set to the upper limits derived from the FUSE analy- 
sis (again without oxygen) and also obtained a good fit (Fig. 12). 
The same result is obtained from a pure H model. 


and they are given in number ratios relative to H. We also com- 
pare line profiles from diffusion models with observations. 

Carbon and oxygen. From the absence of Chi A 1175 A, we 
derive an upper limit of C/H = 1 .0 x 10 -7 . The stratified models 
overpredict the line strengths (Fig. 9). Note that all computed 
line profiles have zero rotation velocity. 

Except for the Ovi AA 1031.91, 1037.61 A resonance dou- 
blet, no O line is identified. Because of a possible ISM contri- 
bution, only an upper limit of O/H = 1 x 10 -6 can be derived. 
The stratified models fit the observation best for logp = 8.2. A 
smaller gravity results in too strong features. 

Silicon, phosphorus and sulfur. From the absence of Siiv 
AA 1066.63, 1 122.48, 1 128.33 A, P v AA 1 1 1 7.98, 1 128.01 A, and 
Siv AA 1062.66, 1072.96, 1073.508 A, upper limits of Si/H = 
5.0 x 10“ 9 , P/H = 1.0 x lO” 10 and S/H = 1.75 x 10“ 8 are de- 
rived (Figs. 10, 1 1). In all cases, the diffusion models predict too 
strong line profiles. 

Iron and nickel. No lines are identified corresponding to up- 
per limits of Fe/H < 1.0 x 1(T 6 and Ni/H < 1.0 x KT 6 . 


5.4. Chandra 

No Fe or Ni lines are detected in the Chandra spectrum, giving 
upper abundance limits of Fe/H and Ni/H < 1.25 x 10 -7 derived 
from homogeneous Fe+H and Ni+H models. These limits are 
almost one order of magnitude more stringent than the FUSE- 
derived results. The observed overall SED can be reproduced 
successfully with both homogeneous and stratified models. The 
interstellar H and He column densities were kept fixed to the 
EUVE-derived values. For the stratified models, however, there 
is still enough Fe and Ni levitated that lines should be detectable 
in the Chandra spectrum, even at a higher gravity of log g =8.5. 

For further analysis, homogeneous and stratified models 
were calculated that incorporated the same metals as the EUVE 
analyses in the preceding section. Individual spectral features in 
the observation are not recognizable, which serves as an addi- 
tional constraint for the upper limits of the metal abundances de- 
termined in the FUSE analysis. Both kinds of models match the 
Chandra spectrum well. The models displayed in Fig. 13 include 
oxygen, the homogeneous model with an abundance according 
to the upper limit given in Table 4. The vast majority of the weak 
lines in these models stems from O. In essence, a pure H model 
gives an equally good fit. 
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Table 3. T & % and log g of GD 246 as determined by previous analyses and the present work (see also Fig. 14). 


No. 

Teff [K1 

logg [cm/s 2 ] 

Method and Model Atmospheres 

Authors 

1 

60100 

7.72 

LTE 

Balmer 

H 

Vennes et al. (1997) 

2 

58 700 

7.81 

LTE 

Balmer 

H 

Finley et al. (1997) 

3 

59000 

7.80 

LTE 

EUVE 

H+metals 

Wolff etal. (1998) 

4 

53 100 

7.85 

NLTE 

Balmer 

H 

Napiwotzki et al. (1999) 

5 

54000 

7.80 

LTE 

Lyman 

H 

Wolff etal. (2001) 

6 

56000 

8.20 

NLTE 

EUVE 

H+metals stratified 

Schuh et al. (2002) 

7 

52402 

7.89 

NLTE 

Lyman 

H+metals 

Barstow et al. (2003a) 

8 

51300 

7.91 

NLTE 

Balmer 

H+metals 

Barstow et al. (2003a) 

9 

54400 

7.90 

NLTE 

Balmer 

H 

Liebert et al. (2005) 

10 

57007 

7.82 

LTE 

Balmer 

H 

Koester et al. (2009) 

11 

54150 

8.00 

NLTE 

Balmer 

H 

Limoges & Bergeron (2010) 

12 

56160 

7.98 

NLTE 

Balmer 

H 

Gianninas et al. (2011) 

13 

57 000 

7.60 

NLTE 

Lyman 

H 

this work 

14 

55 000 

7.90 

NLTE 

Chandra 

H+metals stratified 

this work 
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Fig. 13. Best-fit models (thick lines) for the Chandra spectrum of 
LB 1919 (thin lines). Top : homogeneous model, bottom: diffusion 
model. Both models are computed with H, C, O, Si, P, and S. The ho- 
mogeneous model has abundances equal to the upper limits given in 
Table 4. 

5.5. Mass and distance 

To determine the mass from T e g and loggr, we use DA evolu- 
tionary tracks by Althaus & Benvenuto (1998). We chose the 
tracks with metallicity z = 10 -3 and hydrogen-envelope frac- 
tional mass of M h /M* = 10 -4 and derive M = 0.66+j} p A/ 0 . 
Taking the flux calibration of Heber et al. (1984), we com- 
pute the distance, using the extinction-corrected visual magni- 
tude that is determined from the observed magnitude V = 16.41 
to m Vo = V - 2.174 ■ E(B - V ) with E(B - V) = 0.030. 
With Teff = 56000 ± 2000 K, log g = 7.9 ± 0.3, and H v = 
7.33 x 10 -4 erg cm -2 s -1 Hz -1 (model Eddington flux at 5454 A) 
we find (see, e.g., Rauch et al. 1994) 

d = 7.1 1 x 10 4 ^Hy M io°- 4m ''o- 1 °S9 = 321 pc. (4) 

5.6. Summary on LB1919 

Our fit to the Lyman lines and the Chandra spectrum gave 
T e ff = 56000 ± 2000 K, loggr = 7.9 ± 0.3. The temperature 
is significantly lower (~10 000K) than that derived from previ- 
ous Balmer-line fits. This result is surprising because Barstow 
et al. (2003a) had established that Lyman-line analyses for DAs 


T ett / kK 

60 57 55 52 



log T,„ / K 

Fig. 14. Results from previous analyses of GD 246 (Table 3). Results 
obtained by LTE analyses (open circles) are compared to NLTE anal- 
yses (filled circles) and the analysis with stratified NLTE model atmo- 
spheres (filled square). The triangles mark the two results of this paper, 
obtained by fitting the Lyman lines (full symbol) and the Chandra SED 
(open symbol). 


exceeding T e ff ~ 50 000 K yield systematically higher tempera- 
tures than Balmer-line analyses, i.e., the opposite that we found 
for LB 1919. In their sample, however, there is also one excep- 
tion to this general trend. Very similar to LB 1919, PG 1342+444 
has T eff = 66 750 K and 54 308 K derived from Balmer and 
Lyman lines, respectively (see also Barstow et al. 2002). 

Metal lines are not detected in the FUSE and Chandra 
spectra. We derived upper abundance limits with homoge- 
neous models. The respective models fit the SEDs observed by 
Chandra and EUVE. Stratified models, on the other hand, over- 
predict metal abundances. Reducing the abundances in these 
models to bring them into agreement with observations re- 
quires an excessively high gravity (logg > 8.5), which is 
not compatible with the Balmer and Lyman lines. In essence, 
this confirms the result obtained from the EUVE analysis by 
Landenberger-Schuh (2005). All observations are compatible 
with a pure H atmosphere. 

6. GD246 

6.1. Previous investigations 

Earlier determinations of T et f and log g of GD 246 
(WD 2309+ 105) are summarized in Table 3 and Fig. 14. 
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Table 4. Element abundances in LB 1919 and GD 246 as determined from UV spectra with homogeneous models with 56 000 ± 2000 K and log g = 
7.9 ± 0.3 for LB 1919 and 57 000 ± 2000 K and log g =7.6 ± 0.3 for GD 246. 


WD 

Element 

Abundance (X/H) 

Uncertainty 

Ions 

Diff. 

Literature 

LB 1919 

C 

<1.0 x 10“ 7 

5.0 x 10“ 8 

Cm 

> 



O 

<1.0 x 10“ 6 

2.0 x 10 -6 

Ovi 

> 



Si 

<5.0 x 10“ 9 

5.0 x 10“ 9 

Sim+iv 

> 



P 

<1.0 x lO" 10 

1.5 x lO" 10 

Pv 

> 



S 

<1.75 x 10“ 8 

7.5 x 10“ 9 

S rv+vi 

> 



Fe 

<1.0 x 10“ 6 

1.0 x 10 -6 

Fe v 

> 



Ni 

<1.0 x 10“ 6 

1.0 x 10“ 6 

Niv 

> 


GD 246 

C 

<3.2 x lO" 8 

1.5 x 10“ 8 

C m 

> 



O 

<6.0 x 10“ 8 

5.0 x 10“ 8 

Ovi 

> 

1.6 x 10 -7 “ 


Si 

1.2 x 10“ 7 

2.0 x 10“ 8 

Si iv 


5.0 x 10 -8i 

3.2 x 10- 8c 

1.2 x 10 _7 ° 


P 

3.0 x 10“ 9 

2.0 x 10“ 9 

Pv 

> 

6.3 x 10- 9c 
7.5 x 10- 9i 


S 

5.0 x 10“ 9 

2.0 x 10“ 9 

S vi 

> 

<3.0 x lO" 76 


Ge 

5.0 x 10“ 9 

1.0 x 10“ 9 

Geiv 

< 

2.5 x lO" 9 ' 7 


Fe 

<1.0 x 10“ 6 

1.0 x 10“ 6 

Fe v 

> 

<2.0 x lO" 56 






> 

<1.2 x 10- 7 ° 


Ni 

<1.0 x 10“ 6 

1.0 x 10“ 6 

Niv 

> 

<1.3 x 10- 7 “ 


Notes. Identified ions are noted in the fifth column. Column six indicates the tendency of the abundance of a diffusion model to be stronger (>) or 
weaker (<) than the observed line strength. Literature values for the abundances are given in the last column. 

References. (o) Barstow et al. (2003b). (i>) Wolff et al. (2001). (c) Chayer et al. (2001). M) Vennes et al. (2005). 


The most recent study by Gianninas et al. (2011) is probably 
the most reliable NLTE Balmer-line analysis because it utilizes 
new Stark broadening data (Tremblay & Bergeron 2009). (We 
use the same data in our analysis.) All analyses are based 
on Balmer- and/or Lyman-line profiles with the exception of 
two, which are based on the EUVE SED. One of them (Schuh 
et al. 2002) uses stratified NLTE models, and it results in an 
exceptionally large gravity. 

A small number of lines of trace metals could be identified in 
the FUSE and HST spectra of GD 246. Beside elements known 
to appear in WDs like Si, C, and P, Vennes et al. (2005) identified 
for the first time lines of germanium (Ge iv) in a HST observa- 
tion and determined log (Ge/H) = -8.6 ± 0.2. 

6.2. FUSE and HST 

Pure hydrogen models were used to fit the Lyman lines. We 
found T eff = 57 000 ± 2000 K and log# = 7.6 ± 0.3 (Fig. 15). 
Also plotted are the line profiles of one of our models that has 
parameters that are essentially equal to those found by Barstow 
et al. (2003a) from their Lyman-line analysis (Teir = 52 000 K, 
log# = 7.9). It can be seen that the differences between the two 
sets of profiles are rather subtle and the fit quality is similar. We 
verified for some model parameters that the influence of metal- 
line blanketing on the line profiles is negligible. 

We then determined metal abundances. Some of them can be 
compared to the results by other authors (Table 4). 

Carbon. No C line is detectable. From the absence of the 
C in A 1 175 A multiplet, we find an upper limit of C/H = 3.2 x 
10 -8 . All stratified models predict strong C lines, indicating that 
they clearly overpredict the C abundance. 

Oxygen. The result for O is similar. We derive an upper limit 
of O/H = 6 x 10 -8 from the absence of a photospheric O vi res- 
onance doublet. The stratified models strongly overpredict O. 



Fig. 15. Fit to the Lyman lines of GD 246 (thin line). Overplotted are 
pure H models with T^s = 57 000 K, log g = 7.6 (thick line), and with 
Teff = 52 000 K, log g = 7.9 (dashed line). 



Fig. 16. Fit to P v and Si iv lines of GD 246 with a homogeneous (thick 
line) and a diffusion model (dashed line). The abundances in the ho- 
mogeneous model are Si/H = 1.2 x 10 -7 , P/H = 3 x 10 -9 , and S/H = 
5 x 10“ 9 . 

Silicon. Si iv lines are seen at AA 1122.48, 1128.33 A (Fig. 16). 
We find Si/H = 1 ,2(±0.2) x 10 -7 , that is identical to the result of 
Barstow et al. (2003a). For the diffusion models, the lines fit at 
Teff = 55 000 K and log # =7.9. 
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Fig. 17. Fit to Geiv A 1 189.03 A in GD 246 (thin line). Overplotted are 
a homogeneous model (thick line, Ge/H = 5 x 10 -9 ) and a stratified 
model (dashed line). 

Phosphorus. Fitting the Pv AA 1117.98, 1128.01 A resonance 
doublet (Fig. 16) gives P/FI = 4(±2)x 10 -9 . The stratified models 
predict too strong lines. 

Sulfur. From the absence of the S vi AA 933.38, 944.52 A reso- 
nance doublet, we find S/H < 5(±2)x 10 -9 . The stratified models 
produce excessively strong lines (Fig. 16). 

Iron and nickel. Fe and Ni lines are not detected in the UV 
range. We derive upper limits of 1 x 10 -6 for both species. Fe 
and Ni are overpredicted by the stratified models. 

Germanium. Our fit to the Geiv AAl 188.99, 1229.81 A reso- 
nance doublet gives Ge/H = 5(±1) x 10 -9 (Fig. 17). This is 
about twice the value derived by Vennes et al. (2005) for a 
56 000 K model. The difference can be traced back to NLTE 
effects (Sect. 3.2.1). The diffusion models drastically underpre- 
dict the line strengths. At T e g = 57 000 K and log <7 = 7.6 no 
Geiv lines are visible. 

6.3. EUVE 

Two types of models were calculated. Both include the metals 
identified in the FUSE observation (the homogeneous models 
with the determined abundances, and C and O set to their up- 
per limits). For one type of models, Fe is added, and the other 
type includes Fe plus Ni (setting Fe/H = 1.25 x 10 -6 and Ni/H = 
7.9x 10 -8 in the homogeneous models). Figure 1 8 shows that the 
stratified models with Fe as heaviest element provide a better fit 
to the observation than the models also including Ni. A grav- 
ity of log g = 8.2 is needed, which is significantly higher than 
the result of the Fyman-line fits. The homogeneous models fail 
completely in reproducing the EUVE SED, whis is in agreement 
with the results of Landenberger-Schuh (2005). 

6.4. Chandra 

Vennes & Dupuis (2002) detected Fe lines in the Chandra spec- 
trum of GD 246. They are identified in Fig. 19. Their NLTE anal- 
ysis indicated Fe/H = 3 x 10 -7 . 

Figure 20 displays our best fits using the homogeneous and 
stratified model types as employed for the EUVE analysis. In 
general, as in the EUVE case, models excluding Ni fit better. The 
final homogeneous model fits worse than the Ni-truncated strati- 
fied model, and with an unrealistically low gravity of log g = 7.3. 
The best-fit (Ni-truncated) stratified model gives T e fF and log g 
that is in accordance with the most recent NLTE B aimer-line 
analysis of Gianninas et al. (2011); see Fig. 14. From the ho- 
mogeneous models, we confirm the Fe abundance derived by 
Vennes & Dupuis (2002). 




Fig. 18. Best-fit models to the EUVE spectrum of GD 246 (thin line). 
Top: stratified models; bottom: homogeneous models. Solid lines: mod- 
els including Fe; dashed lines: models include Fe plus Ni. 

6.5. Mass and distance 

With T eff = 57 000 ± 2000 K, log g = 7.6 ± 0.3, H v = 7.52 x 
10 -4 erg cm -2 s _I Hz -1 , m v = 13.09, E(B - V) = 0.0011, we 
derive M = 0.54 / 947 M e and d = 94^ pc. 

6.6. Summary on GD246 

Our Lyman-line fit gave T e s = 57 000 K and log g = 7.6. The 
Balmer-line fit by Gianninas et al. (2011) resulted in a simi- 
lar temperature (56 160K) but at a significantly higher gravity 
(7.98). With our parameters, the stratified model compared to 
the FUSE data overpredicts C, O, Fe, Ni, P, and S, but under- 
predicts Ge. That result remains essentially valid for stratified 
models with increased gravity (log g = 7.9). 

No acceptable fit to the EUVE spectrum was achieved. 
The gravity of the best fit homogeneous model is too low 
(log g = 7.3), while that of the best-fit stratified model is too high 
(log g = 8.2). It must be stressed that in both cases the formally 
best fits are very poor. 

For fitting the Chandra spectrum, the stratified models are 
only better than the homogeneous models when Ni is arbitrar- 
ily removed. The best-fit stratified model has then log 17 = 7.9, 
which is closer to the Balmer-line than to the Lyman-line result. 
However, the fits deteriorate when Ni is included. But even when 
Ni is excluded, the best-fit stratified model is not satisfactory. We 
conclude that stratified models do not predict correctly the indi- 
vidual element abundances. This conclusion is corroborated by 
the FUSE analysis. 

7. PG 1520+525 

7.1. Previous investigations 

The first analysis of PG 1520+525 (WD 1520+525) based 
on optical spectra was performed by Werner et al. (1991) 
and gave T e s = 140000 ± 20 000 K and log g = 7.0 ± 1.0. 
A first indication that the temperature exceeds 140 000 K 
was derived from the shape of an EUVE spectrum (Werner 
et al. 1996). Utilizing UV spectra, Dreizler & Heber (1998) 
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Fig. 19. Chandra spectrum of GD246 (thin 
line; bin size 0. 1 A) and model spectrum 
(shifted upward by 150 counts; Fe POS lines 
only). 



Fig. 20. Best-fit models (thick lines) to the Chandra spectrum of 
GD246. Top : homogeneous model with abundances according to 
Table 4. Bottom : stratified model. Both models have no Ni. 
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Fig. 21. Three models with different T e g fitted to the Chandra spectrum 
of PG 1520+525. Model abundances are given in Table 5, log g = 7.5, 
loglV(Hi) = 20.0. 


Table 5. Atmospheric parameters of PG 1520+525 and PG 1159-035. 



PG 1520+525 

PG 1159-035 

Teff [K] 

150000+ 10000 

140000 + 5 000 

log 9 (cgs) 

7.5 + 0.5 

7.0 ±0.5 

He 

0.43 

0.33 

C 

0.38 

0.48 

O 

0.17 

0.17 

Ne 

0.02 

0.02 

Mg 

0.006 


Fe 

0.0013 

0.0013 


Notes. PG 1520+525: T e g, log g from Dreizler & Heber (1998) and our 
work. Abundances (in mass fractions) from Werner et al. (2004b), ex- 
cept for Mg (this work) andFe (Werner et al. 2011). PG 1159-035: from 
Werner et al. (201 1). 

found T e ff = 150000± 10000K and log g = 7.5 ±0.5. These pa- 
rameters, as well as the He, C, O, and Ne abundances (Werner 
et al. 2004b), are given in Table 5. 

7.2. Chandra analysis 

Most counts in the Chandra spectrum are located in 
the 100- 120 A region. At longer wavelengths, the ISM absorbs 
the photospheric flux, while at shorter wavelengths, the flux de- 
creases due to photospheric O vi bf absorptions. Although S/N 


is rather poor, it is possible to identify individual lines from O vi 
and Ne vn (Fig. 22). The identification of Mg lines is ambiguous. 
For a detailed list of possible lines in this wavelength range we 
refer to the analysis of H 1504+65 (Werner et al. 2004a). 

We have calculated a small grid of models for T e ff = 
140 000 K, 150 000 K, 160 000 K, and log g = 6.4, 7.0, 7.5. We 
kept fixed the He/C/O/Ne abundance ratio but varied the Mg 
abundance from solar to ten times solar in order to derive an 
upper limit. 

The fit is hardly affected by the particular choice of log g. 
Instead, Teff is the most sensitive parameter. Mainly judging 
from the spectral slope, we find the best fit at T e ff = 150000 ± 
10 000 K adopting log g = 7.5 ± 0.5 (Fig. 22). This confirms the 
previous temperature determination from UV/optical results, but 
the Chandra spectrum allows no further reduction of the error 
bar. 

A change in the Mg abundance from solar to ten times solar 
affects the strength of the lines only marginally, and it is not clear 
which abundance is more likely to fit the observation (Fig. 23). 
The flux of the models is also affected in regions where no Mg 
lines are expected. This can be explained by the occurrence of 
autoionization features in the bf cross sections. These resonances 
can be strong and narrow, looking like a common absorption line 
(Werner et al. 2004a). The wavelengths of these resonances are 
uncertain and could thus be a possible explanation for unidenti- 
fied features in the observed spectrum. 
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Fig. 22. Line identifications in PG 1520+525. Overplotted are two mod- 
els with different T e t fi 150 000 K (thick line), 160 000 K (dashed line). 
Model abundances are given in Table 5, log g = 7.5, logA(Hi) = 20.0. 



Fig. 23. Detail of the Chandra observation of PG 1520+525 (thin line), 
overplotted with two models with solar (dashed line) and ten times solar 
(thick line) Mg abundance (r e ff = 150 000 K, log g = 7.5). 

7.3. Summary on PG 1520+525 

Our analysis of the Chandra spectrum gives an independent con- 
firmation of the result on T e (j derived from previous UV spec- 
troscopy (150000± 10000 K, assuming loggr = 7.5 ± 0.5). 

As outlined in Sect. 2.2, we may use this result to draw con- 
clusions on the location of the blue edge of the GW Vir instabil- 
ity region, which is confined by the non-pulsator PG 1520+525 
and the pulsator PG 1159-035. The position of both stars in 
the Teir -log g diagram is displayed in Fig. 24. Also shown there 
is the theoretical location of the blue edge predicted by Quirion 
et al. (2007) assuming He/C/O mass fractions of 0.4/0. 4/0. 2 in 
the stellar envelope, which is within error limits identical to the 
abundances found in both stars (Table 5). Besides these element 
abundances, the exact location of the edge depends on the rnetal- 
licity; therefore, it is shown here for z = 0.0 and z = 0.007. 
The edge moves to higher effective temperatures with increas- 
ing z. The z = 0.007 line is appropriate for the two stars con- 
sidered here, because we had found solar Fe abundance in them 
(Werner et al. 2011). It can be seen that the theoretical position 
of the edge runs between the location of both stars, very close to 
PG 1159-035. 

Also shown in Fig. 24 are other PG 1 159 stars in which pul- 
sations were looked for (taken from the compilation in Werner 
& Herwig 2006). The location of nonpulsators redward of the 
indicated blue edge is possibly due to their different chemical 
composition compared to PG 1 159-035, because the composi- 
tion affects the position of the instability strip (Quirion et al. 
2007). 

8. Summary and conclusions 

We analyzed the Chandra spectra of three hot WDs. Two of them 
are H-rich WDs and the other one is a H-deficient PG 1159 star. 

8.1. The DA white dwarfs LB 1919 and GD246 

The primary aim of our Chandra observation of LB 1919 was to 
explain its relatively low metallicity compared to similar objects 
as suggested by previous analyses of EUVE spectra. It turned out 
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Fig. 24. Pulsating (filled circles) and nonpulsating (empty circles) 
PG 1159 stars in the logT e ff - logg diagram. Evolutionary tracks are 
labeled with the respective stellar masses in M e (Miller Bertolami & 
Althaus 2006). The red edge (short dashed line) of the instability region 
(Quirion et al. 2006) and two blue edges (long dashed lines, Quirion 
et al. 2007) are shown. The upper and lower blue edges are for z = 0 
and 0.007, respectively. 

that no metal features were detected in the Chandra spectrum. 
The same result was found from our analysis of a FUSE spec- 
trum. In essence, all data are compatible with the assumption 
that LB 1919 has a pure hydrogen atmosphere. This is in conflict 
with our stratified models. In them, the vertical run of individual 
element abundances is computed from the assumption of equi- 
librium between gravitational downward pull and radiative up- 
ward acceleration, i.e., the metal abundances are not free param- 
eters but computed as functions of T e ff and log g. They predict 
that significant amounts of light and heavy metals should be ac- 
cumulated and readily detectable in the atmosphere of LB 1919. 

A few other hot DAs with similar parameters that also do 
not have detected metals are known (Barstow et al. 2003b); 
one famous example is HZ43A (T e ff = 51 000 K, log g = 7.9, 
Beuermann et al. 2006). The reason for the purity of their atmo- 
spheres is completely unknown. It may be speculated that these 
peculiar DAs have no heavy-element reservoir that is assumed to 
be present in the radiative levitation calculations. Barstow et al. 
(2003b) consider whether the possible depletion of those reser- 
voirs by selective mass-loss could be responsible for that phe- 
nomenon or whether the progenitors were metal-poor. They ar- 
gue that, for radiation-driven winds, mass loss should be lower 
in the pure H stars than in those containing heavy elements, be- 
cause the wind is driven by metal lines. Furthermore, these WDs 
are local disc objects, so it seems unlikely that any of the pro- 
genitors could have been metal poor. 

We have also analyzed an archival Chandra observation of 
the hot DA GD246. Our work was motivated by the fact that 
this is the only DA that exhibits individual metal-line features 
in the soft X-ray range. As for LB 1919, we utilized chemically 
homogeneous as well as stratified models. While we in principle 
expect that the stratified models are a better representation, the 
Chandra observation cannot be fitted satisfactorily by either type 
of model. A better fit by stratified models is only achieved when 
nickel is removed artificially. Otherwise, the atmospheric opac- 
ity becomes much too large because of a strong overprediction 
of nickel. In comparison with the FUSE spectrum, most met- 
als are overpredicted, but at the same time, another one (Ge) 
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is underpredicted. In the latter case, one can speculate that the 
model atoms are still too small, not having enough line transi- 
tions to absorb photon momentum. The overprediction of most 
metals, however, is more of a problem. We can only speculate 
that an additional physical mechanism, ignored by our models, 
is affecting the equilibrium abundances, e.g., (selective) mass- 
loss or weak magnetic fields. 

The complete failure of our model atmospheres, homoge- 
neous or stratified, to fit the EUVE observation reinforces the 
conclusion that some physics is missing. The same conclusion 
was drawn by Barstow et al. (2003b) from the fact that UV spec- 
tra from DAs with similar T e ff and log g display a wide variety 
of metal abundances. These authors considered the possibility 
that accretion of interstellar or circumstellar matter could cause 
this variety. Recent investigations have revealed the frequent in- 
cidence of dust debris disks around WDs (e.g., Farihi 2011) that 
are composed of disrupted planetary material. In order to test 
the influence of accretion of such material on the observed pho- 
tospheric abundances, this effect would have to be included in 
the diffusion/levitation models in a manner performed (e.g., by 
Koester 2009) in cooler WDs in which radiative levitation is neg- 
ligible. However, it is difficult to understand how an external 
supply of material would resolve the problem of present diffu- 
sion models generally overpredicting metal abundances. 

Our analysis of the two DAs comprised the derivation of ef- 
fective temperature and gravity from the Lyman lines, and the 
resulting parameters were compared to published results from 
Balmer-line analyses. We confirm previous findings that con- 
flicting results are obtained from the UV and optical line pro- 
file analyses. In particular, Balmer-line temperatures are often 
significantly lower than Lyman-line temperatures (Barstow et al. 
2003a). We found that LB 1919 is another rare case where the 
opposite effect was found. While this obviously points to a short- 
coming in the WD atmosphere models, there is no indication as 
to which physical ingredient is treated inadequately or is miss- 
ing. These uncertainties in the derived atmospheric parameters 
add to the problems in a detailed quantitative comparison, ele- 
ment by element, of observed metal abundances with predictions 
from radiative levitation models. 

8.2. The PG 1 159 star PG 1520+525 

The Chandra spectrum of the nonpulsator PG 1520+525 was 
used to constrain its effective temperature and to compare it with 
that of the pulsator PG 1159-035. They both confine empiri- 
cally the blue edge of the GW Vir instability region for a par- 
ticular envelope composition. The position of the edge predicted 
by the nonadiabatic pulsation models of Quirion et al. (2007) 
is consistent with the spectroscopic results. This is strong proof 
for the predictive power of these models. The interior structure 
of PG 1159 stars that is inferred by their usage appears to be 
very reliable, so that corresponding asteroseismologic analyses 
are based on solid ground. 
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